Strain-stiffening in random packings of entangled granular chains 
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Random packings of granular chains are presented as a model polymer system to investigate the 
contribution of entanglements to strain-stiffening in the absence of Brownian motion. The chain 
packings are sheared in triaxial compression experiments. For short chain lengths, these packings 
yield when the shear stress exceeds a the scale of the confining pressure, similar to packings of 
spherical particles. In contrast, packings of chains which are long enough to form loops exhibit 
strain-stiffening, in which the effective stiffness of the material increases with strain, similar to 
many polymer materials. The latter packings can sustain stresses orders-of-magnitude greater than 
the confining pressure, and do not yield until the chain links break. X-ray tomography measurements 
reveal that the strain-stiffening packings contain system-spanning clusters of entangled chains. 

PACS numbers: 83.80.Fg, 81.70.Tx, 62.20.mm, 61.41. +e 



Most materials, ranging from crystalline metals to piles 
of granular material, become weaker the further they are 
strained. On the other hand, many polymeric materials 
are known to strain-stiffen, in which the effective stiffness 
of the material increases as the material is strained fur- 
ther pQ. Theories suggest strain-stiffening could depend 
on many factors including chain stiffness, density, tem- 
perature, strain rate, and in particular on structures such 
as entanglements between different chains [2H5]. How- 
ever, it has not been possible to directly measure en- 
tanglements in experiments because in polymers these 
structures occur on very small scales. In this letter we 
present a new experimental approach to investigate the 
role of entanglements in strain-stiffening using a model 
system of granular chains consisting of millimeter-scale 
beads connected by flexible links. Such a macroscopic 
system has advantages over molecular polymer systems 
for investigating entanglement. First, the macroscopic 
size allows for imaging to measure the precise positions 
of each particle and link in the structure. Second, we can 
isolate entanglement effects from temperature and strain 
rate dependent effects because the macroscopic chains 
have no inherent time scales due to Brownian motion or 
relaxation. 

Granular chains have been considered as model poly- 
mers in previous studies to characterize the packing 
structure around the jamming transition for chain pack- 
ings [6-8 . Such chains form tight loops when packed 
which defines a characteristic loop size and persistence 
length in analogy to polymers. The free volume in the 
packing was found to increase with and level off when 
the chain length exceeded the persistence length, analo- 
gous to the behavior of the glass transition temperature 
for polymers [7 j. Similarly, tumbling strings have been 
considered as a macroscopic model polymer system, in 
which knotting was found to exhibit statistics analogous 
to thermodynamic systems [9J. While these works char- 
acterized the structures of macroscopic model polymer 
systems, the stress response of such systems have not 
previously been characterized. In this letter, we report 




FIG. 1: (a) An individual chain with length N — 20 beads. 
The minimum loop size of £ = 8 beads is highlighted by the 
red circle, (b) Cylindrical packing of chains with N ^> £ = 8 
inside a thin latex membrane used for triaxial compression 
measurements. Many tight loops near the minimum loop size 
can be seen in the packing through the membrane. 



stress measurements to demonstrate a more direct con- 
nection with polymer materials; specifically that they ex- 
hibit similar strain-stiffening. We use x-ray tomography 
to measure the precise packing structure and identify en- 
tanglements to demonstrate a quantitative connections 
between entanglements and strain-stiffening. 

As a model material, we used macroscopic granular 
chains consisting of hollow spherical brass beads shown 
in Fig. [TJi. The beads are flexibly connected into chains 
by enclosing dog-bone-shaped brass links. These connec- 
tions have essentially zero stiffness for small bend angles, 
but they have a maximum bend angle beyond which the 
chains will not bend without plastic deformation. This 
defines a minimum loop circumference £ the chains can 
bend into, as shown in Fig.[l^L. We used two sets of chains 
with beads of slightly different diameters of 1.9 mm and 
2.1 mm which causes them to have different minimum 
loop circumferences of £ = 8 and £ = 11 beads, respec- 
tively. In contrast to polymer chains, these chains have 
a non-linear response to bending as described above, and 
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FIG. 2: Stress r vs. strain 7 for packings of chains with length 
N w 10 4 . The shaded bands represent 1 standard deviation 
based on 5 repeated measurements. Strong strain-stiffening is 
seen as the positive curvature in r(j). Inset: r(j) for N = 1 
(unlinked beads). 



they are far too large to experience Brownian motion. 

We measured the mechanical response of granular 
chain packings using triaxial compression [10]. The 
grains were poured into a flexible elastic membrane in 
a cylindrical shape with solid end caps on the top and 
bottom, as shown in Fig.fl}}. An Instron materials tester 
was used to quasi-statically push the top end cap down- 
ward while measuring its height and the required force 
for samples with diameter and initial height of 50 mm. 
This average compressive stress r is calculated as the 
force over the cross-sectional area of an end cap and the 
compressive strain 7 is calculated as the height change 
(positive downward) over the initial height. 

We show in Fig. [2] the stress-strain relation t(j) for 
a packing of chains with length N ~ 10 4 beads, much 
longer than their minimum loop circumference £ = 8. 
Enormous strain-stiffening can be seen as the region of 
positive curvature of t(j). A striking feature is that 
the stress reaches the order of 10 MPa before the pack- 
ing fails, which is indicated by the flattening of the 
stress-strain relation where the packing is able able to 
shear without supporting additional stress. This ulti- 
mate strength is on the order of 1/10 the strength per 
density of the solid brass that makes up the chains (this 
packing is about 20% brass by volume) and comparable 
to polymer rubbers. Such high strength per density in a 
pourable, shapable granular material could have unique 
engineering applications. In contrast, unconnected gran- 
ular materials typically fail when the shear stress exceeds 
the confining pressure at the boundary in the absence of 
interparticle attractions [10]. In our experiments, the to- 
tal confining stress coming from gravity and the elastic 
membrane is only about 10 kPa, which is the maximum 
stress scale reached for N — 1 (unlinked beads) in the 
inset of Fig. [2] (the positive slope for 7 > 0.05 matches 
the contribution of the membrane stiffness as it is de- 
formed). The chain packings exceed this stress by a fac- 
tor of ~ 10 3 . The strength of polymer materials is partly 
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FIG. 3: Stress r vs. strain 7 for different chain lengths N 
shown in the key. (a) For chains of minimum loop circum- 
ference £ = 8. (b) For £ = 11. In both panels, the data for 
N > £ (cool colors) group into a band which exhibits strong 
strain-stiffening. 

attributed to the extensional or bending strength of in- 
dividual chains. Indeed, during the compression of the 
long chains in Fig. [2] we started to hear chains break at 
the point where the maximum stress was reached, and 
after each measurement was done we counted approxi- 
mately 10% of the links in the chains to be broken. This 
breaking at such stresses suggests that - similar to poly- 
mer systems - the material strength of the chains them- 
selves contributes and the packing can only shear and fail 
by breaking links in the chains, which would require the 
packing structure to produce some self-confinement. 

To investigate the transition from granular behavior 
in the limit of unlinked beads with N = 1 and strain- 
stiffening in the limit of large A 7 ", we show stress-strain 
curves for packings of chains with different lengths N in 
Fig. [3] Chains of minimum loop circumferences £ = 8 and 
£ = 11 are shown in panels a and b, respectively. It is seen 
in both panels that the data fall into two distinct bands: 
the upper band of data with chain lengths A" > £ exhibits 
strong strain-stiffening, while the lower band of data with 
A" < £ exhibits typical granular strain-softening or much 
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minimum loop circumference £ 


8.0 ±0.5 


11.0 ±0.3 


min. length for strong strain-stiffening 


8.0 ±0.5 


11.0 ±0.5 


min. length for majority clusters 


6.5 ±1.0 
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FIG. 4: (a) Examples of entanglements for chain lengths 
N — 5 (left) and N — 20 (right). The green chains are 
entangling the purple chains. The shaded region indicates 
the entanglement manifold connecting the ends of the chains 
(m = N — 1). (b) Average number of entanglements per con- 
nector as a function of connector separation m for £ = 8. The 
key indicates the chain lengths N . (c) Maximum size of en- 
tangled clusters as a fraction of the total number of chains 
in the sample. The dotted line indicates the onset of strong 
strain-stiffening, which is only found to occur when the con- 
nector separation is large enough to that there are entangled 
clusters consisting of the majority of chains in the system. 
The inset image shows an example of a cluster of 4 chains 
with 5 entanglements between them for N — 10. 



weaker strain-stiffening. The fact that the transition 
length between these two bands scales with and is about 
equal to the minimum loop size £ in both cases suggests 
a relationship between strain-stiffening and whether the 
chains form loops in the packings. 

To investigate the role of entanglement in this relation- 
ship, we used x-ray tomography to precisely reconstruct 
the positions of each bead and link in three dimensions. 



TABLE I: Comparison of different length scales. The chain 
length required for strong strain- stiffening is just above that 
required for most of the chains to be entangled in a single 
cluster, and consistent with the minimum loop circumference 
of the chains. 



For these experiments, we used aluminum chains with 
bead diameter 2.5 mm and £ = 7.5 in cylindrical pack- 
ings with diameter and height of 35 mm each. We show 
results for chain lengths N = 5, 10 and 20 each com- 
pressed to a strain of 7 = 0.2. 

We define an entanglement to occur if one chain wraps 
partially around another. Motivated by string knotting 
analysis [9], we first draw an imaginary line called the 
connector between any two beads of the entangling chain 
separated by a distance of m beads along the contour 
of the chain. We then define an entanglement mani- 
fold as the minimal two-dimensional manifold bounded 
by the connector and the links that connect the inter- 
vening beads. If another chain crosses through this man- 
ifold, it is counted as entangled by the entangling chain. 
Examples of entanglements along with the manifolds are 
shown in Fig. [4^i. These particular examples would not 
be counted as entanglements based on typical polymer 
algorithms where entanglements are counted if chains 
catch on each other if they are contracted [11 . In con- 
trast with elastic polymers, our granular chains are highly 
constrained from significant rearrangements due to their 
high packing density, so we expect these examples are still 
able to produce significant constraints against shearing 
from entanglement. Regardless, the majority of our en- 
tanglements would satisfy the typical requirements, and 
it is usually expected that different algorithms for mea- 
suring entanglement tend to produce qualitatively similar 
statistics [TT] . 

We count all of the entanglements for each possible 
connector between beads and plot the average number of 
entanglements per connector as a function of separation 
distance m in Fig. [4)3 [12]. Regardless of chain length 
TV, for connector separations of m < 3, no entanglements 
are observed because the size of the beads restricts other 
chains from fitting inside even the largest possible en- 
tanglement manifolds. For larger m, the average num- 
ber of entanglements per connector increases monotoni- 
cally with m. Since chains that are longer compared to 
their minimum loop size can be bent further, the area 
of the entanglement manifold tends to similarly increase 
with m and it is more likely in a random packing that 
other chains will cross this manifold to become entan- 
gled. Thus, the entanglement manifold method provides 
an intuitive way to understand why the probability of 
entanglements increases with separation or chain length. 

The increase in the number of entanglements per con- 
nector with separation distance seen in Fig. [4b could 
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help explain why only the longer chains strain-stiffen. 
The longer chains with N = 10 and 20 which strain- 
stiffen reach more than 1 entanglement per connector 
for larger connector separations, while the shorter chains 
with N = 5 which do not strain-stiffen have no more than 
0.04 entanglements per connector. Since the probability 
of entanglements as a function of separation distance is 
similar for different iV, it seems likely that the reason 
that the N = 5 chains do not entangle more is simply 
that they are too short compared to the minimum loop 
size £ = 8 for them to bend much and thus have entan- 
glement manifolds of large enough size to result in many 
entanglements. 

If entanglements are responsible for constraints that 
prevent failure by shear, then at minimum the entangle- 
ments would have to connect into system-spanning clus- 
ters to avoid weak points in the packing that could still 
shear. We identify a cluster as subset of chains in entan- 
glements that are transitively connected. An example of 
a cluster of four entangled chains is shown in the inset of 
Fig. [4J3. The largest cluster size found for each N and m 
is plotted in Fig. For the strain-stiffening chains with 
N = 10 and 20, the largest cluster size jumps sharply to 
over 50% of the chains in the packing at m — 6 and 7, 
respectively. This is just below the chain length TV = 8 
required for strong strain-stiffening (these length scales 
are summarized in Table [i]). These cluster sizes also large 
enough that they can span the system in several direc- 
tions. This is necessary for there to be enough constraints 
from entanglements in those clusters to fully constrain 
the packing against shear. In contrast, the N = 5 chains 
which do not strain-stiffen also do not have large enough 
clusters of entanglements to span the system. 

To summarize, we demonstrated that random pack- 
ings of granular chains can be used as a model sys- 
tem that exhibits strain-stiffening as in polymer mate- 
rials, with the additional ability to precisely measure the 



structure of entanglements. These chains exhibit strain- 
stiffening if they are longer than the minimum loop cir- 
cumference they make in packings, which allows them to 
form system-spanning clusters of entanglements. Since 
the strength of granular materials is usually limited by 
the scale of the confining pressure at the boundaries, 
we propose that for the chain packings to reach the 
orders-of-magnitude higher strength, the confinement at 
the boundary must be functionally replaced by a self- 
confinement due to the system-spanning clusters of en- 
tanglements. Strain-stiffening then results if these en- 
tanglements tighten up under strain, eventually locking 
up to prevent further shear unless the links break. This 
picture contrasts with recent theoretical arguments that 
strain-stiffening in polymers can be explained without 
entanglement. However, those arguments only explained 
much weaker strain-stiffening in which the stress increase 
is less than an order of magnitude [H [5], so no addi- 
tional confining stress mechanism was necessarily in those 
cases. The questions of whether entanglement is the only 
way to get strong strain-stiffening, and how to quantita- 
tively predict strain-stiffening from entangled structures, 
remain open for future work. 
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